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AINS LAND SURVEYOR SYSTEM WITH 
REPROCESSING, AINS-LSSRP 

This non-provisional application claims priority from U.S. provisional 
application serial number 60/252,862 filed November 22, 2000. 

5 L Field of the Invention 

The subject invention relates to inertial navigation, and more particularly to 
aided inertial navigation systems (AINS) based Land Surveyor Systems which use an 
AINS as a navigational reference and which make it possible to survey forested areas 
where GPS signals may be missing for long intervals of time or indefinitely due to 

10 foliage and or a dense tree canopy. Land surveyor systems that are not aided with a 
radio positioning system such as GPS as a result of signal blockage due to a tree 
canopy tend to accumulate position error linearly as a function of time. 

2. Background of the Invention 

AINS technology originated in the late 1960's and found application on 

15 military navigation systems. An example or the many books on the subject is the text 
by George Siouris, titled "Aerospace Avionics Systems, A Modern Synthesis", 
published by Academic Press, published in 1993. 

Traditional methods of surveying used laser theodolites, which required access 
to lines of sight between positions to be surveyed. A line of sight is not always 

20 available in forested areas. A line of sight was typically obtained in the past by 

bulldozing corridors 3-4 meters wide through the forest along the lines or paths to be 
surveyed. Governments discourage this wasteful and environmentally destructive 
practice by imposing stumpage fees on destroyed trees and other forms of penalties for 
environmental damage. The advent of a precise GPS has provided an alternative 

25 method of land surveying; however, a Land Survey System based exclusively on a 
precise GPS receiver is limited to areas where the sky is visible. Some tropical rain 
forests have canopies that are so dense that precise GPS cannot be used at all. 
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AINS Land Surveyor 

An AINS land surveyor does not require access to the sky and can be operated 

under a dense tree canopy. A single surveyor using an INS land surveyor can map a 

forest or jungle area by walking a path among the trees thereby avoiding the need to 
5 cut trees to establish a survey lane as would be required using theodolites or a precise 

GPS. Use of an AINS land surveyor therefore results in reduced cost, and the 

environmental impact of the survey is substantially reduced. 
^ To reduce drift and obtain acceptable accuracy, the surveyor using an AINS- 

vP based land surveyor systems brings the INS to a complete rest about every 1-2 minutes 

y§ 10 for a period of 5-30 seconds. This is called a zero-velocity update (ZUPD). A 

Kalman filter uses each ZUPD to zero the INS velocity error and partially calibrate out 
£ inertial sensor errors. The position error drift with periodic ZUPD's is on the order of 

, 0.5-2 meters per kilometer, depending on the quality of the inertial sensors and on the 

Q frequency of ZUPD' s. The requirement for ZUPD's is often an inconvenience, since 

\f; 15 it limits the surveyor's production. 

13 The Kalman filter in an AINS based Land Surveyor System is also coupled to 

receive measurement inputs from one or more external and independent sources which 
are processed to further refine or aid the navigation solution. Aiding information can 
include aiding signals from sources such as a precise GPS receiver, a Doppler Radar, 

20 or from an odometer or a precise pedometer. In the survey of an area permitting an 
unobstructed view of the sky, a precise GPS provides the simplest aiding signal. 
Several different accuracy levels are associated with GPS performance and signals. 
C/A GPS implies uncorrected GPS, and provides 10-20 meter 3D position accuracy. 
A GPS can be coupled via a radio modem link to a base GPS receiver, which has a 

25 precisely known position. The base GPS measures the errors in the signals being 
received and forwards corrections via the model to the precise GPS aboard the 
navigator. An RTCM-corrected differential GPS uses industry standard RTCM 
differential corrections from a dedicated base receiver or a differential corrections 
service such as the U.S. Coast Guard or Omnistar to obtain 0.5-1 .0 meters position 

30 accuracy. Real-time kinematic (RTK) GPS uses differential GPS data from a 
dedicated base receiver to obtain 0.05-0. 1 meters position accuracy. 
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The invention system allows a surveyor to carry the invention AINS-based 
Land Surveyor System with Reprocessing along a predetermined path or line and to 
locate and record the position of stakes that the surveyor positions in accordance with 
a pre-planned grid of locations. The stakes identify the locations of sensors and 
5 explosive charges that are used in a seismic survey. The predetermined path begins 
and ends at position fix (PF) locations that are precisely known. The invention AINS- 
based Land Surveyor System with Reprocessing then uses the known PF locations at 
O each end of the path and a smoother algorithm to reprocess the stake locations. 

y§ Seismic Exploration 

jLJ? 10 Seismic exploration has as its object, the production of a multi-dimensional 

=P map of the geological structure over an area below the ground for the purpose of 

identifying valuable oil, gas and mineral deposits. A seismic survey uses acoustic 
[7 interferometry to perform the multi-dimensional subterranean mapping. A 

N geophysicist provides a pre-plot map with grid locations of the desired positions of the 

□ 15 noise sources and the geophones over the space to be explored. The noise sources are 

multiple phased dynamite explosions on a 2-dimensional grid pattern. The sound 
waves from the charges are reflected by the different geological strata, and received by 
an array of geophones on a separate 2-dimensional grid connected to recording 
devices. 

20 The multi-dimensional geological map is generated by post-processing the 

recorded data. The typical error specification for the position of a noise source or a 
geophone is one meter horizontal and 0.5 meters vertical. 

A backpack-borne AINS based system is made and sold by the assignee of this 
application, Applanix Corporation, at 85 Leek Cresent, Richmond Hill, Ontario, 

25 Canada L4B 3B3. The system is called the Position and Orientation System for Land 
Survey (POS/LS). In operation, a surveyor walks a survey path or trajectory carrying 
the POS/LS as a backpack. Such survey trajectories often pass through areas where 
GPS signals are not available. The POS/LS navigates though such GPS outage areas 
in a dead-reckoning mode with as little position drift as possible. Typically the 
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surveyor moves from one known position to another, and "ties-down" or "fixes" the 

POS/LS position at the known positions. 

The path followed by a surveyor is typically a zigzag pattern of parallel 

seismic lines, each 1000-5000 meters long with the stakes positioned every 50-200 
5 meters on the lines. A control survey is sometimes performed to verify the accuracy 

of the seismic survey. A control survey is formed by a number of short survey paths 

between known positions with traverse legs being approximately perpendicular to the 
pi principal seismic or grid lines. Intersections of the control paths with the seismic lines 

^ form the control points. The control point position accuracies should be significantly 

# 1 0 better than the seismic line position accuracies. 

!f SUMMARY OF THE INVENTION 

This invention teaches the use of an optimal smoothing algorithm and a 
M method of reprocessing the data accumulated and stored from an AINS-Land Surveyor 

y ; System with Reprocessing (AINS-LSSRP) between known position fixes or "tie- 

Jr 15 downs". 

A first embodiment of the AINS-LSSRP comprises an AINS that provides a 
sequence time-indexed present position values in response to the LSSRP being moved 
from a first known position value or PF at the start of a survey interval to a second 
known position value at the end of the survey interval. The LSSRP has a Position and 
20 Orientation (POS) Computer Subsystem (PCS) coupled to receive and store the 

sequence of time-indexed present position values as the surveyor moves from the first 
known position to the second known position. The PCS has a reprocessing computer 
and program means for processing the indexed present position values with a 
smoothing algorithm to provide indexed and adjusted present position values for at 
25 least some of the indexed present position values recorded as the LSSRP was moved 
from and between the first known position value at the start of the survey interval and 
the second known position value at the end of the survey interval. 

In a more particular alternative embodiment, of the LSSRP, the AINS uses a 
Kalman filter responsive to at least two sources of aiding signals, and the PCS has an 
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aiding signal selector algorithm characterized to select the most accurate aiding signal 
for use by the Kalman filter from all available aiding signals. 

In yet another more particular embodiment of the LSSRP, the AINS uses a 
reprocessing computer and program means smoothing algorithm that is used is a 
5 Modified Bryson-Frazier smoother (MBFS). 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram of an aided INS system; 

Fig. 2 is a block diagram of the Aided Inertia! Navigation System - Land 
Surveyor System with Reprocessing (AINS-LSSRP); 
10 Fig. 3 is a schematic characterization showing a top graph of real-time position 

error growing linearly with time; a lower graph shows reprocessed position solutions 
with a maximum average error at a location approximately between position fixes; 

Fig. 4a is a graph of real-time position error* the system using ZUPD's each 
minute to aid the AINS between position fixes; 
15 Fig. 4b is a graph of real-time position error, the system using ZUPD's each 

minute to aid the AINS between position fixes, the positions being reprocessed; 

Fig. 5 is a flow chart for a computer program for storing time position values 
as the system moves between position fixes and for reprocessing the real-time position 
values at the end of a position fix; 
20 Fig. 6 is a flow chart for signaling the start of recording smoother data at the 

beginning of a survey interval and for subsequently closing the smoother data file at 
the end of the survey interval; and 

Fig. 7 is a flow chart for signaling the start of running the smoother filter at the 
end of the previous survey interval. 

25 DETAILED DESCRIPTION OF THE INVENTION 

The invention taught herein improves the accuracy of reported stake positions 
by combining the benefits of an AINS, such as the AINS used in an Applanix Position 
and Orientation System for Land Survey (POS/LS) with the steps of onboard data 
reprocessing. The combination forms an AINS- Land Surveyor System with 
30 Reprocessing (AINS-LSSRP) represented by block 42 in Figure 2. In operation and 
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use, a surveyor begins a survey interval or leg from a position fix (PF) location having 
a precisely known location. The surveyor initializes the present position of his AINS- 
LSSRP using differential GPS or a pre-surveyed marker. The surveyor then follows 
predetermined survey lines or paths that run for several kilometers in accordance with 
5 desired grid locations on a map. The surveyor may follow a zigzag pattern where the 
survey lines go into tree covered areas and where GPS contact is sometimes lost In 
some rain forests, such as those in Indonesia, GPS can be lost for an entire day. The 
m surveyor uses the AINS-LSSRP to position stakes at the predetermined grid locations 

W and on reaching a second PF, the surveyor resets or "ties down" the AINS-LSSRP 

ytj 10 position with the precise position coordinates of the PF. The AINS-LSSRP then 

launches the reprocessing function, which reprocesses and adjusts data gathered in the 
+5 course of the traverse between the two PF locations including the recorded positions 

s ' of the stakes using a smoothing algorithm. The adjusted positions of the stakes are 

[7 reported as a product of the survey. 

^ 15 Figure 1 shows a generic Aided INS architecture. The elements within 

3 phantom box 20 represent an inertial navigation system (INS), having an inertial 

measurement unit (IMU) 22 and an inertial navigator 24. The inertial navigator is 
typically mechanized using a digital computer and software for processing signals 
from and for delivering signals to the IMU 22. 
20 The IMU 22 comprises a triad of accelerometers (not shown) that measure 

total specific force to restore a proof mass, the force being proportional to 
acceleration, and a triad of gyros (not shown) that measure total angular rate. The 
IMU also provides the supporting process and interface electronics (also not shown) 
that convert and output the inertial acceleration and angular rate signal data in a digital 
25 format. The inertial navigator is a mechanization of Newton's equations of motion, in 
software, typically running on the aforementioned digital computer. To do this, the 
INS 20 typically aligns or transforms the signal data from its package or vehicle 
navigation coordinate frame into a fixed and earth referenced coordinate system, such 
as a North, East and Down directions. In some mechanization's, the horizontal axes is 
30 aligned to a heading other than North/East and the heading offset angle is called the 
wander angle. A typical ground alignment requires the INS to be stationary for 5-15 
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minutes so that the INS can measure the gravity vector with its accelerometers to 
establish Down alignment and with a knowledge of the latitude, measure the 
horizontal component of the earth rotational rate vector with its gyros to establish the 
North direction and or the wander angle. 
5 An aided inertial navigation system (ADSfS) 26 typically comprises an INS 20, 

one or more aiding sensors such as the GPS receivers 28, a Doppler Radar 30 or a 
distance measuring instrument (DMI). A Kalman filter 34 and a mechanism such as 
the error controller 35 process corrections to the INS 20. 
Notation 

10 The following notation is used in this description of the invention. Lower case 

letters denote vectors. For example, x denotes a vector. Upper case letters denote 

matrices. For example, H denotes a matrix. 

Subscripts on vectors and matrices are used to indicate a discrete sequence 

number that correspond to a sample time. For example, the matrix Hk denotes that the 
15 time of validity of the matrix H is sample time ft, and hence is the &-th sample in a 

discrete time series of//. 

"E" denotes the expectation operator. See A. Papoulis, "Probability, Random 

Variables and Stochastic Processes", McGraw-Hill, 1965 or any other reference on 

random variables for its definition. 

20 Kalman Filter 

The Kalman filter 34 is a recursive minimum-variance estimation algorithm 
that computes an estimate of a state vector based on constructed measurements. The 
measurements typically comprise computed differences between the inertial 
navigation solution elements and corresponding data elements from the aiding 

25 sensors. For example, an inertial-GPS position measurement comprises the 

differences in the latitudes, longitudes and altitudes respectively computed by the 
inertial navigator and a GPS receiver. The true positions cancel in the differences, so 
that the differences in the position errors remain. A Kalman filter designed for 
integration in an INS with aiding sensors estimates the errors in the INS outputs and 

30 aiding sensors. INS errors typically include: inertial North, East and Down position 
errors; inertial North; East and Down velocity errors; inertial platform misalignment 
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in the local level vector and wander angle errors; accelerometer biases and gyro 
biases. Aiding sensor errors can include the following: GPS North, East and Down 
position errors; GPS carrier phase ambiguities and DMI scale factor error. 

The Kalman filter therefore performs as an optimal linear minimum variance 
5 estimator that operates on a time series of measurements to estimate the state vector of 
a linear stochastic multivariable model of the form: 

where 

x k is the n-dimensioned state vector at discrete time 
10 <£(A+1 ; k) is called the transition matrix, or the nxn system dynamics state 

transition matrix. It is the discrete version of the system matrix <X> and is analogous to 
the F matrix that appears in the text by G. Siouris referenced above. 

Jl is an n-dimensioned vector of process noises with variance Qk 

where 

is E[j2 k julj] = Q k S } (2) 

The subscripts k and k+j are sample counts or indices. The samples are valid 
at discrete times t of k) and t of (k+j). The dependency on time t is dropped and only 
the sample indices of k and j are retained to identify the discrete samples. 

The jX symbol is a way of specifying white noise. Equation (2) above 

20 represents the expected value of a noise vector at one time multiplied by the noise 
vector transpose at a subsequent time, the product being equal to the covariance 
matrix Q times the delta function which basically says that the expected value is non- 
zero only when the two times are the same. In other words, noise vectors sampled at 
different times are completely uncorrelated. The arrow above //* indicates that it is an 

25 /i-dimensional vector of process noises with variance. The process noise is not a 
known error. Process noise is the driving noise in the stochastic model for the 
dynamics that the Kalman filter is designed to estimate. The state vector is therefore a 
vector of random variables whose randomness requires the use of a Kalman filter. 
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If a process has no noise, then it should be possible to predict every fixture 
value of that process based on some initial value. If the first value is known and the 
system dynamics state transition matrix is O (Phi) that describes how that process 
proceeds, then in the absence of noise, every future value can be predicted exactly. A 
Kalman filter is not necessary in such systems. 

z k is the m-dimensioned measurement vector at discrete time &. In operation, 

it is equal to a position fix at discrete time 

Hk is the mxn measurement model matrix, or measurement matrix. 
fjk is the m-dimensioned vector of measurement noises with variance: 



As above, terms that begin with "E[ ]" imply an expected or average value, 
where S k is the discrete delta function given by: 
8k = 1 when k = 0, Sk = 0 otherwise. 

<1>(A;-H; k) describes the transition of the dynamic system from sample k at 
time t k to sample k+l at time t k +j. A text by T. Kailath, "Linear Systems", Prentis- 
Hall, Inc. 1 9 1 00 provides a further discussion. Such a model is often used to model 
the behavior of a noisy multivariable process. A Kalman filter is used to estimate the 
internal dynamic parameters of the process based on noisy measurements of the 
process. This is the conventional method of modeling and estimating the errors of an 
INS system such as that shown in Figure 1. The Kalman filter is constrained to causal 
operation, hence its estimate of a state vector is the optimal estimate based on all past 
and current measurements. The following equations describe the standard Kalman 



Extrapolation or time update: 

Equations 4 and 5 below extrapolate the state vector x k from sample time t^j 
to sample time f*. The state vector x k is updated at sample time t k using 
measurements valid at time k and again for the next sample time 4+/. 




(3) 



filter: 



(4) 



P; =<S>(k:k-l)P;_&{k:k-\) r +Q, 



(5) 
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Equations 4 and 5 show that the transition matrix® transitions the state 
estimate of the state vector x k from sample k-1 to sample k. The transition matrix 
&(fc, k-l) describes the system dynamics in a discrete setting. 

Measurement update: 

**=4~+**f* (6) 
£ k =z t -H k x- k (7) 
P;={l-K k H k )P k - (8) 
K k =P k -H T k S- k l (9) 
S k =H k P k 'H T k +R k (10) 
where xl is the priori (before measurement update) estimate of the state vector, 
xl is the posteriori (after measurement update) estimate of the state 

vector, 

£ k is the Equations 4 and 5 innovations vector. The term "innovations" is 
an established and accepted term in Kalman filtering that describes the 
difference between the measurement vector %k and the model for the 
measurement vector before the update is processed, 
is the priori estimation error variance-covariance (VCV) matrix, 
is the posteriori estimation VCV matrix, 

is the Kalman gain matrix, 
is the innovations VCV matrix. 
The VCV matrices /£* , P£ and Sk are propagated in parallel with the state 

vector x£*~ using equations (5) and (8). The discrete system dynamics are defined per 

equation (1). The VCV matrix is then propagated. 

The estimation error VCV matrix describes the Kalman filter's predicted 
estimation error propagation statistics as follows: 

p;'-*E[(r k '--x k )(x;'--x k ) T ] (ii) 



S k 
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The diagonal elements of equation (1 1) are the predicted estimation error 
variances for the individual state vector elements. If the model of equation (1) exactly 
describes the system dynamics and measurements, then the estimation error VCV 
matrix exactly describes the estimation error statistics of equation (11). The model of 
equation (1) is at best an approximation for a nonlinear system; therefore the 
estimation error VCV matrix is interpreted for each case. 

The minus and plus exponents in the above equations imply values 
respectively occurring before and after the Kalman filter measurement update at a 
given time index k. A (-) superscript indicate a before measurement update. A (+) 
superscript indicates an after measurement update. A (+/-) indicates that the equation 
holds for either case. This notation is used in a related text titled "Applied Optimal 
Estimation", MIT Press 1994 (12th printing 19 1 12) by A. Gelb, (editor). 

The theoretical value of the covariance matrix is the expected value of the 
error vector times the error vector transpose. The product of two parameters is the 
cross correlation or the autocorrelation between those two values. If they are 
uncorrelated from one sample to the next then the autocorrelation is zero except at the 
time when the two samples are valid or at a point where the two samples are taken at 
the same time. They are equal to the variance if it is a scalar or to the covariance if it 
is a vector* 

At the conclusion of each AINS system navigational solution, a time blended 
position is the result. If the Kalman filter is designed to interface to a smoother, it 
writes the following variables to a file which the smoother later reads: 
<3>(&-H ; k) transition matrix from k to A+l 
Hk measurement matrix at iteration k 

Kk Kalman gain at iteration k 

S/c Kalman innovations covariance at iteration k 

innovations vector at iteration L See page 31 7 in Gelb's 
"Applied Optimal Estimation" text referenced below. 

Kalman estimated error state at iteration k 

Kalman generated estimation error VCV matrix at iteration k 
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To implement a Kalman filter, the designer first chooses or codes the 
parameters <£(ZH-1 ; k), Hk and Rk so that the Kalman filter properly models the 
process being estimated. If the process being estimated is the INS error process, then 
the parameters k) and Qk are derived from one of several known INS error 

5 models. Error models are discussed in the above referenced text by George Siouris, 
titled "Aerospace Avionics Systems, A Modern Synthesis", published by Academic 
Press, published in 1993, The designer then designs the measurements z k that the 
Kalman filter will process, and from these measurements, the filter computes and 
implements the measurement model parameters Hk and For a GPS-aided INS, the 
10 measurements z k are typically differences between the INS and GPS position and 
velocity data. Hk describes the measurement geometry with respect to the state vector 
x k , and Rk describes the measurement noise variance or relative weighting of each 

scalar measurement in the measurement vector z k . The steps in Kalman filter design 

are discussed in several texts, one of which is mentioned above, i.e., the text by A. 
15 Gelb titled "Applied Optimal Estimation", MIT Press 1994 (12th printing 19 1 12) by, 
(editor). 

Equations 4 through 10 above describe the full mechanization of the Kalman 
filter, given the above model parameters and measurements. There are also several 
references on numerically stable and efficient implementations of these equations. A 
20 popular implementation is that disclosed as the Bierman implementation described in 
the text titled "Factorization Methods For Discrete Sequential Estimation" by G. J. 
Bierman and published by Academic Press in 1977* 

Optimal Smoothing Algorithms and Methods 

A smoothing algorithm is defined here to be any algorithm that computes an 
25 estimate of a quantity or collection of quantities assembled into a vector at a given 
estimation time using measured data having times of validity before and after the 
estimation time. 

An optimal fixed interval smoother computes the optimal estimate based on all 
past, current and some future measurements. It is characterized in the text by A. Gelb, 
30 referenced above, as a variance weighted-weighted combination of a Kalman filter 
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running forward in time on a data time series and a Kalman filter running backward in 
time on the same data time series. The fixed interval smoother operates 
on all available measurements to compute the optimal estimate at each sample time 
based on all measurements during a defined process interval. 

Different smoother formulations can be obtained by the algebraic manipulation 
of the basic forward-backward Kalman filter smoother. The following equations 
describe one such formulation that is developed further in the text title "Applied 
Optimal Control" by A. E. Bryson Jr. and Y. C. Ho, published by Wiley and Sons, 
1975. The smoother is referred to as the modified Bryson-Frazier smoother (MBFS). 

Backwards extrapolation: 

& x =Q> T {k:k-\)X; (12) 

A U = ^ r ( k : * -1) A;a>(* :k~l) (13) 
Adjoint measurement update: 

l" k ={l-H k K k )l\ -H k S~ k l Z k (14) 
A; =(l-H k K k f hl(l-H k K k )+H T k S k H k (15) 

where: 

X k ~ is the a priori adjoint state vector, 
A* is the a priori adjoint VCV matrix, 
X£ is the a posteriori state vector, and 
A£ is the a posteriori adjoint VCV matrix. 

The input data to the MBFS at the A-th iteration of the MBFS comprises the 
following data items that a Kalman filter designed for interface to a MBFS computes 
internally and writes to data files: 

<t>(k; k-1) transition matrix from iteration k-1 to iteration k 
Hk measurement matrix at iteration k 

Kk Kalman gain at iteration k 

Sk Kalman innovations covariance at iteration k 

% k innovations vector at iteration k. Seepage 317 in Gelb*s 

"Applied Optimal Estimation" text referenced above. 
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x£ Kalman estimated error state at iteration k 

Kalman generated estimation error VCV matrix at iteration k 
The smoothed state and estimation error VCV matrix are obtained as follows: 

xt = r k -p;x; = % -p;X t + (16) 
5 p;=^/-a^)=p;(/-a;/> + ) (17) 

Note that smoothing with the a priori or a posteriori states and VCV matrices 
yields the same smoothed state and VCV matrix. This is the mathematical expression 
of the fact that the smoothed estimate at any estimation time is the optimal estimate 
based on all measurements, before and after the estimation time. 

10 Application of Optimal Smoothing to the AINS-LSSRP 

A smoother designed for an AINS-LSSRP comprises two components: the 
smoothed estimation operation and the position solution correction operation. 

The smoothed estimation operation uses the MBFS described previously that 
is made specific to the AINS component of the AINS-LSSRP primarily through the 
15 design of the stochastic state and measurement models in Equation 1 . The Kalman 
filter is designed to estimate the errors in the INS and siding sensors. This a well- 
understood practice among navigation practitioners, and is described in G. Siouris 
referenced above and reviewed here briefly. 

The elements of the state vector x h in Equation 1 comprise the INS navigation 

20 errors (position, velocity and attitude errors), inertial sensor errors (accelerometer and 
gyro biases and possibly other errors) and aiding sensor errors. The AINS-LSSRP 
requires an improved position solution only, hence the state vector is given the 
following representation: 
Sr 



^remaining 



(18) 



25 where Sr is the AENS position error 

^remaining ls ^ e sub-vector of error states other than Sr 
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The elements of the measurement vector z k in Equation 1 comprise 

differences between the positions and velocities computed by the INS and 
corresponding positions and velocities from the aiding sensors* The Kalman filter 
thus compares components of the INS navigation solution with corresponding 
components from the aiding sensors, and thereby estimates the errors in the INS 
navigation solution and the aiding sensors. 

The AINS then corrects the INS navigation algorithm with the estimated INS 
navigation errors. The closed-loop between the INS and the AINS Kalman filter 
shown in Figure 1 achieves regulation of the INS errors to consistency with the aiding 
sensor errors. If the adding sensors include a GPS receiver, then the AINS navigation 
errors will be consistent with the GPS navigation errors. In particular, AINS error 
regulation suppresses the INS position error drift, and the Kalman filter estimates and 
thereby calibrates the inertial sensor errors such as the accelerometer and gyro biases. 
By using a Kalman filter to estimate the INS errors, the AINS achieves a theoretically 
optimal (in the least-squares sense) navigation error based on all past and current 
aiding sensor data. 

The MBFS described in Equations 12 to 15 operates on data generated by the 
AINS Kalman filter and written to data files, and computes a smoothed error state and 
VCV described in Equations 16 and 17 that have improved accuracy over the 
estimated error state from the AINS Kalman filter. The improvement results from the 
fact that the MBFS computes the optimal estimate (in the least-squares sense) at an 
iteration k based on all aiding sensor data before and after the iteration k. The MBFS 
thus computes estimates of the INS position, velocity and attitude errors, called the 
smoothed INS errors, that the AINS error regulation mechanism was unable to correct 
because of the causal constraint on access to aiding sensor data. The smoothed error 
state is partitioned as follows compatible with Equation 18. 




6f s 



(19) 



'remaining 



where Sr s is the smoothed estimate of AINS position error 



'remaining 



is the sub-vector of smoothed error states other than 5f s 
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The position solution correction operation corrects the AINS position solution 
by simple subtraction of the smoothed AINS position errors from the AINS position 
solution. The actual mechanization depends on the AINS position coordinate format, 
which can be one selected from the following possible formats: geographic (latitude, 
longitude, altitude), Universal Transverse Mercatur (UTM) (Northing, Easting, height) 
or earth-fixed Cartesian (X,Y,Z components of position vector The following 

example corrects an earth-fixed Cartesian position vector r£ to obtain the corrected 

or reprocessed position vector ff ♦ 



(20) 



10 Zero-Velocity Update (ZUPD) 

The AINS-LSSRP represented by block 42 in Figure 2 is brought to a 
complete rest periodically, typically every 1-2 minutes, for a period of 5-30 seconds. 
Each rest period is called a zero-velocity update (ZUPD). A Kalman integration filter 
uses these zero velocity observations to zero the INS velocity error and to null inertial 

15 sensor rate errors. An AINS-LSSRP (Land Surveyor System with Reprocessing) can 
reduce position error drift with periodic ZUPD's to a value of 0.5-2 meters per 
kilometer depending on the quality of the inertial sensors and the time between 
ZUPD's. However, the requirement for ZUPD's limits the surveyor's production. 
Referring to Figure 3, survey lines are called traverses. The surveyor starts 

20 each traverse interval from a known position, referred to in Figure 3 as a PF interval, 
an indexed position fix at a mapped survey grid location. The surveyor then follows 
the survey lines while positioning stakes at required map grid locations. The surveyor 
continues along the traverse and tries to end up at another known position and does a 
"tie down" where he establishes the next PF or position fix. "Tying down" means that 

25 the surveyor arrives at the second known position and enters the known coordinates of 
the position fix into this navigation system. The position error that has accumulated 
over the traverse is reset to zero 
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The reprocessing feature of the invention activates the processor's smoother 
algorithm at each "tie-down". Reprocessing enables the system to use the recently 
stored data to calculate the errors and to correct the position of stakes located in the 
course of the traverse interval between position fixes. Having corrected the locations 
5 of the stakes for a completed traverse, the system is free to retain or delete the 
previously stored time-indexed position values if memory constraints require the 
storage space used for the previously stored time-indexed position values be made 
available for use in the course of the next survey interval. 

The time indexed solutions are position values that the AINS computes and 

10 logs in a mass memory unit as the surveyor is traversing a survey line between 
position fixes. As a tie down location is reached, the operator pushes a button or 
otherwise initiates the reprocessing function. In the alternative, the operator could 
transfer the known position values for the tie down as a pre-recorded waypoint, push 
the button and keep going. The computer can navigate with time navigation on a 

15 second traverse while running the smoother on the recorded positions from a first 
traverse. 

As the smoothing algorithm completes the smoothing process, it outputs a 
smoothed error table and surveyed positions with improved accuracy for the data 
gathered in the previous traverse. The smoothed error table is a close estimate of the 
20 position error for each time position in the previous traverse. Curve 58 in Figure 3 
schematically depicts the re-processed position error from a completed survey 
interval. 

Figure 2 is a block diagram of an embodiment of the invention ABSfS-LSSRP 
42 that uses a generic smoothing algorithm to re-process the aided inertial navigation 

25 data. The AINS-LSSRP 42 comprises an AINS 26 with aiding signal sources therein 
as shown in Figure 1, a processing computer called the POS computer subsystem 
(PCS) 44 comprising a STORAGE DEVICE 46, a reprocessing computer and 
program means 45, and an aiding signal selector 43. Input/output interface device is 
coupled to the PCS 44 and allows the survey to input data into the AINS 26, 

30 reprocessing computer and program means 45 and aiding signal selector 43 as 
required. 
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Input/output interface device 48 is typically a keyboard and/or a display, with 
which the operator uses to initiate reprocessing at a known tie-down or position fix. 
The interface device has a mode control bus 36 coupled to the aiding signal selector 
input 51, a position fix entry bus 37 coupled to the AINS 26 from output 39 and an 
5 input/output bus 38 coupled from output port 40 to the reprocessing computer and 
program means input 41 . The mode control bus 36 and the input/output bus 38 
provide input and output control signals and position values and smoothed position 
values as required between the respective functional blocks for control, display and 
yj data transfer purposes. 

10 The AINS 26 provides a series of real-time position values via bus 27 to the 

W reprocessing computer and program means 45 for transfer and storage in arrays of 

s p memory locations in mass storage memory 46 via data bus 49. The time-indexed 

J" position values are output during each traverse between "tie-downs'* via bus 27 to the 

^ reprocessing computer and program means block 45 and via data bus 49 to mass 

15 storage memory 46. 

2"l The reprocessing computer and program means 45 reprocesses and filters each 

N * set of time-indexed position values. The aiding signal selector 43 is programmed or 

controlled by an input from the surveyor via the input/output interface device 48 via 
bus 36 to direct the AINS 26 to select the most accurate of the aiding signals 
20 available, such as a precise GPS aiding position signal. 

Input/output interface device 48 receives position values (PF values) from the 
surveyor and transfers the values to the AINS 26 and to the reprocessing computer and 
program means 45. The reprocessing computer and program means 45 is 
programmed with the algorithm that mechanizes a smoother filter such as the filter 
25 characterized by equations 12 through 17 above. 

The aiding signal selector 43 is typically a computer program, a sequential 
machine or software program, or an operator-controlled software switch that selects 
the highest accuracy time-indexed aiding position values available and couples the 
selected aiding signal values to the AINS 26 or directs the AINS 26 as to which signal 
30 source to select. By way of example, if the radio modem drops out, the GPS position 
accuracy can drop from 0.05 meters to 20 meters. If the AINS-LSSRP has maintained 
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a position accuracy of better than a few meters, the aiding signal selector would direct 
the AINS to treat a radio modem link loss as a complete GPS outage. Alternatively 
the operator may be suspicious of GPS accuracy and would then manually deselect 
GPS data processing by the AINS. 
5 The AINS-LSSRP will not use C/A quality GPS data because it can achieve 

better position accuracy between two position fixes using ZUPD aiding alone. In this 
mode, the AINS-LSSRP navigates as if GPS data is not available in the absence of a 
^ radio modem link. Re-processing is important in this case even though the receiver 

"B has clear access to GPS C/A quality signals from the GPS satellites. Since the radio 

jqi 10 modem can be lost, without notice, at any time, ZUPD's continue to be necessary for 

15 improved accuracy. 

Ms?- 

Figure 2 shows that aiding signal selector is coupled to AINS 26 via signal 
f path 25. It should be understood that all data might be coupled between the all of the 

M blocks shown via a common bus connecting the functions characterized in Figure 2. 

e 15 The reprocessing computer and program means 45 reads the set of time-indexed 

L f position values and smoother values from the Kalman filter in INS 26 and uses the 

start position fix or tie-down and end position fix or tie-down and all of the time- 
indexed position values from the previous related traverse to calculate an estimated 
error for each time-indexed position value or PF for each traverse. The error is then 
20 applied to coirect each corresponding time-indexed position value to provide a set of 
corrected or adjusted time-indexed position fixes. Each filtered time-indexed position 
fix has an accuracy that is improved over the accuracy of the corresponding time- 
indexed position value fix. 

In a more particular embodiment of the AINS-LSSRP 42, the reprocessing 
25 computer and program means 45 for reprocessing and filtering a set of time-indexed 
position values is programmed to be a smoother such as a Modified Bryson-Frazier 
Smoother (MBFS) that operates on all available measurement values to compute 
optimal estimates of position values based on all measurements during a defined 
processing interval or traverse. The MBFS operates on a set of time-indexed position 
30 values for each traverse using the start PF and the end PF values with the position 
values benefiting from all ZUPD's between the start and end position fixes. 
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Figure 3 schematically shows that the real-time position error, top curve 54 
increases linearly with time and reaches a maximum at peak 56 on the last position 
solution before the AINS Kalman filter processes the next position fix. By contrast, 
curve 58 shows that the reprocessed position solution exhibits a maximum average 
error approximately midway between position fixes. A survey interval therefore 
comprises a PF interval plus a subsequent traverse interval plus the next PF interval. 
A survey interval thus begins with the start of a PF interval and ends with the end of 
the subsequent PF interval. 

Figure 4a shows a plot of real-time radial position error data as a function of 
time in thousands of seconds. Figure 4b shows a plot of the reprocessed position error 
data as a function of time in thousands of seconds. ZUPD's were used to aid the 
AINS between position fixes. The reference position for computing the radial error 
was obtained from a precision GPS position solution with 0.05 meter accuracy. 

A Position Fix 

A "position fix started" condition is determined when the operator directs the 
ABSfS-LSSRP to perform a position fix/ZUPD at a pre-surveyed location. The system 
is stationary at the time; therefore, a position fix is a combination of position fix and a 
ZUPD. The operator enters the position coordinates using the input/output interface 
device 48, typically a handheld screen/keyboard. 

As an alternative, a precise GPS calculated position may be transferred via a 
link such as an infrared or optical, acoustic or electrical link in place of an operator 
manually entered position to avoid errors and reduce the operator workload. As 
presently conceived, a manual action from the operator initiates or triggers the transfer 
of a position fix from the GPS to the system at the location of a tie-down. As an 
alternative embodiment, the GPS position is used in place of the manually entered 
position to avoid operator-induced errors. 

Blocks 52 - 60 of Figures 5 shows the steps performed by the generic AINS 
26, described schematically in Figure 1. The steps begin at ENTER block 52. The 
AINS 26 receives IMU data such as accelerations and rotation rates via block 54 at the 
IMU sampling rate, of typically 50-500 Hz and occasionally as high as 1000 Hz. The 
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inertial navigator algorithm 56, operates on the data, to update each IMU record and to 
provide a time-indexed INS navigation solution. 

The decision 58 to run the Kalman filter can be based on time or on 
availability of aiding values. If, for example, a GPS receiver nominally provides 
5 position information every 1 second, then AINS can run the Kalman filter every 1 
second with the arrival of each batch of new GPS data. Alternatively, if only ZUPD's 
are available, the Kalman filter is run at the completion of each respective ZUPD. 

At the conclusion of each cycle, or indexed computation, the Kalman Filter, 
represented by block 34 in Figure 1, and the RUN THE KALMAN FILTER AND 
10 CORRECT THE INERTIAL NAVIGATOR WITH ESTIMATED ERRORS 

(FIGURE 1) block 60 in Figure 5, updates its estimate of the INS errors. The Kalman 
filter provides its estimate to the error controller block 35, on Figure 1, that is used to 
send navigator correction signals to the Inertial Navigator 24. 

Figure 5 shows the RUN THE DR FUNCTION FOR SURVEY INTERVAL 
15 [k, k+1] (FIGURE 6) block 90, and the RUN THE DR FUNCTION FOR SURVEY 
INTERVAL [k+1, k+2] (FIGURE 6) block 110. Figure 5 also shows the RUN THE 
SEC FUNCTION FOR SURVEY INTERVAL [k, k+1] (FIGURE 7) block 130 and 
the RUN THE SEC FUNCTION FOR SURVEY INTERVAL [k+1, k+2] (FIGURE 7) 
block 150. 

20 The DR and SEC Functions 90, 130 in Figure 5 collectively describe the 

reprocessing function. The DR and SEC Functions run whenever the AINS-LSSRP 
Kalman filter 34 runs. The Kalman filter can process a single position fix 
measurement or a succession of measurements of the same position fix during a 
Position Fix (PF) interval during which the AINS-LSSRP is stationary at a known 

25 position. The Kalman filter will processes ZUPD measurements, precise pedometer 
or other aiding input measurements but not position measurements during a traverse 
interval. The AINS-LSSRP position error will drift at a specified rate. In processing a 
position fix during a PF interval after a traverse interval, the Kalman filter updates its 
estimate of AINS position error and reduces its position uncertainty from the 

30 accumulated position drift during the traverse interval to the uncertainty in the 
position fix, typically one the order of a few centimeters to a few decimeters 
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depending on the accuracy of the surveyed position or of the precise GPS position 
information available. 

Two pairs of blocks 90, 130 and 110, 150 depict two successive PF intervals 
during which the Kalrnan filter processes a succession of measurements of a position 
5 fix. The first time traverse interval is identified as the [k, k+1] interval. The [k, k+1] 
interval starts when the AINS-LSSRP begins to be stationary at a known position and 
the Kalman filter begins to process measurements every n-th second (typically every 
one second) of the known position. The Kalman filter may also process other 
^ measurements such as ZUPD measurements or pedometer measurements if these are 

m 10 enabled as part of the AINS-LSSRP. The [k, k+1] interval ends when the AINS- 

rij LSSRP begins to move and the Kalman filter stops processing position fix 

t; measurements. Interval [k, k+1] can describe a single Kalman filter position fix 

JS measurement, in which case the [k, k+1] interval starts and ends at the same time, or a 

succession of measurements, in which case the interval [k, k+1] starts and ends at 
f7 15 different times. The survey interval [k, k+1] in DR Function 90 and SEC Function 

# 130 Figure 5 extends from the start of interval [k, k+1] at sample "k" to the end of the 

f l [k, k+1] interval at sample "k+1". Likewise, the survey interval [k+1, k+2] in the DR 

Function 1 10 and SEC Function 150 in FIGURE 5 extends from the start of interval 
[k+1, k+2] to the end of interval [k+1, k+2]. The survey intervals [k, k+1] and [k+1, 
20 k+2] will overlap at the end of the [k, k+1] interval. 

Therefore, FIGURE 5 shows two parallel DR and SEC function boxes to 
indicate the possible concurrent execution of two DR and SEC functions for survey 
intervals [k, k+1] and [k+1, k+2]. Each pair of processes, a DR process followed by a 
respective SEC process, is therefore a generic process, each pair being started by a 
25 host process for a given survey interval and then terminated when their operation is 
ended. 

In the alternative, each pair could be the functional object of a generic class 
that is initialized and executed for the survey interval [k, k+1] or [k+1, k+2] assigned 
to them, and then re-initialize for the next survey interval [k+3, k+4] or [k+4, k+5] 
30 assigned to them. 



23 



The PR Function 

FIGURE 6 shows the DR function for the survey interval [k, k+l]. The sub- 
routine of Figure 6 begins when the sub-routine exits block 60 via signal path 86 on 
Figure 5 and enters Figure 6 at START block 92 after which it advances to HAS 
5 INTERVAL [k,k+l] STARTED, decision box 94. If decision box 94 determines that 
the traverse for interval [k, k+l] has started, the routine exits decision box 94 via YES 
path 96 to the HAS INTERVAL [k, k+l] ENDED decision box 98. If the routine 
determines that the survey interval [k, k+l] has not ended, exits via the NO path 100 
3 to the RECORD SMOOTHER DATA box 102. The DR Function routine then 

10 records the time-present position values to a smoother data file in a set of smoother 
JJJ data files that it open on first pass through box 102. The routine continues to record 

4S smoother values exiting box 102 to the DONE box 103 on each subsequent pass until 

f the HAS INTERVAL [k, k+l] ENDED decision box 98 has determined that [k, k+l] 

has ended. On determining that interval [k, k+l] has ended, the routine exits decision 
M; 15 box 98 via the YES path 104 to the CLOSE SMOOTHER DATA FILE box 106 at 

pi which time the routine closes the smoother files in the allocated memory. Each time 

p ~ the routine completes the task of blocks 102 or 106, the routine passes via paths 107 

or 108 respectively to the DONE block 103, after which it follows path 132 in Figure 
5 to the RUN THE SEC FUNCTION FOR SURVEY INTERVAL [k, k+l] (FIGURE 
20 7) block 130. Figure 7 shows the steps in the SEC routine of block 130. The routine 
is entered via signal path 132 to start block 131. Thereafter the DR function of block 
90 for the survey interval [k, k+l] ceases to be active. 

The DR function thus records smoother values from the start to the end of 
interval [k, k+l]. This set of smoother files is then applicable to the survey interval 
25 [k, k+l]. Each successive DR function generates a set of smoother files for each 
successive survey interval. 

Typically survey intervals are defined between sequential start and end 
position fixes (i.e., one after the other), but can be defined to include other position 
fixes between them. The DR function generates files of values to be reprocessed, 
30 containing the navigation values to be corrected and the AINS Kalman filter values 
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required by the smoother algorithm. The DR function generates one file or set of files 
for each survey interval to be reprocessed. 

As the operator enables reprocessing, the DR function writes the values 
required by the smoother to a smoother data file or set of files on the mass storage 

5 memory 46 beginning at a starting position fix interval PF1 and ending at an end 
position fix interval PF2 for a given survey interval. The DR function determines the 
time to start and stop recording when the operator directs the AINS-LSSRP to enter 
positions at PF1 and PF2 via the input/output interface device 48. At PF1, the DR 
function opens the smoother data file(s) and begins recording subject to the AINS- 

10 LSSRP confirming that it is at rest. Data continues to be recorded as the AINS- 
LSSRP moves to the next position fix at PF2 where the operator brings the AINS- 
LSSRP to rest and provides a PF2 signal. The DR function closes the smoother data 
file(s) when the AINS-LSSRP begins to move again. The SEC function then runs the 
smoother on the recorded values and computes smoothed navigation solutions. The 

15 navigation solutions can be for the stake or control positions, or for a timed series of 
values during the survey segment. 

As the AINS-LSSRP continues to move, the DR function continues to write 
values required by the MBFS to the next smoother data file on the mass storage 
memory 46 beginning at a starting position fix PF1 and ending at an end position fix 

20 PF2. The DR function determines the time to start and stop recording when the 

operator directs the AINS-LSSRP to enter positions at PF1 and PF2 via the interface 
device 48. As before, at the end of the PF2 interval, the DR function closes smoother 
values file. 

The SEC Function 

25 FIGURE 7 shows the SEC function for the survey interval [k, k+1]. The sub- 

routine of Figure 7 begins on Figure 5 or Figure 6 when the sub-routine exits block 90 
via signal path 132 and enters Figure 7 at START block 131 after which it advances to 
the HAS INTERVAL [k, k+1] ENDED ? decision box 134. The test for this 
determination can involve a manual signal from the surveyor or an evaluation of 

30 velocity or position signals to determine if the package has again started moving. The 
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SEC Function for the survey interval [k, k+1] will run concurrently with the DR 
Function for the survey interval [k+1, k+2]. The exceptions are the first and last 
position fixes in a survey job. The SEC Function tests for the end of an interval [k, 
k+1] during which position fixes are processed. If the decision box 134 determines 
5 that the interval [k, k+1] is not ended, i.e., the DR function continues to record time- 
indexed position values as the traverse continues, and the SEC Function routine 
advances via NO path 135 to the DONE box 146 and then via signal path 148 to 
ENTER box 52 on Figure 5. 

If the decision box 134 determines that the [k, k+1] interval is ended, the SEC 

10 Function routine advances via YES path 136 to the HAS THE DR FUNCTION 

FINISHED RECORDING AND CLOSED DATA FILE FOR INTERVAL [k, k+1]? 
decision block 138 to determine if the DR Function has recorded a complete set of 
smoother values for the survey interval [k, k+1]. 

If the data set is not complete, the SEC Function routine advances via the NO 

15 path 139 to the DONE box 146, exiting the DONE box via signal path 148 to Figure 5 
to ENTER box 52 for another iteration. 

If the data set exists and is complete, the SEC Function exits decision box 138 
via the YES signal path 140 to the RUN SMOOTHER AND CORRECT POSITION 
SOLUTIONS box 142. Box 142 functions to launch the smoother, which processes 

20 the data set and computes a set of smoothed navigation solutions for the [k, k+1 ] 

survey interval. The smoothing and error correction (SEC) function of box 142 run a 
smoothing algorithm, such as the MBFS of equations 12 - 17 on the smoother data file 
to compute smoothed estimates of position error in the real-time position solutions. 
The smoother then corrects the real-time position solutions to obtain smoothed 

25 position solutions with better accuracy than the real-time position solutions. 

The smoothed navigation solution has the best achievable accuracy with the 
values recorded during the survey interval, and hence is called the best estimate of 
trajectory (BET). When the smoother has completed processing the smoother values, 
the SEC Function routine for the survey interval [k, k+1] exits via signal path 144 to 

30 the DONE box 146. The routine exits the DONE box via signal path 148 to Figure 5 
to ENTER box 52 for another iteration. On exiting box 146, the SEC Function for the 
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traverse of interval [k, k+1] increments the index counter for the next traverse and the 
[k+1, k+2] interval. After exiting block 130, the SEC function for interval [k, k+1] 
ceases to be active. 

As the SEC Function routine completes the SEC computation of the time- 
5 indexed position value data stored during the traverse of interval [k, k+1], all files 
required by the SEC for this survey interval can be deleted to limit the growth of data 
storage space. In principal, this function is not required if data storage space is 
p% unconstrained or otherwise not an issue. 

|f On the next pass, the subroutine operates the same performing the same steps 

yp 10 on Figure 5 from box 52 through box 60. As a result of the index counters being 

Qj incremented from "k" to <4 k+l" and "k+1" to "k+2", the next traverse leaves block 60 

j: via path 88 and starts the next traverse with the RUN THE DR FUNCTION FOR 

SURVEY INTERVAL [k+1, k+2] (FIGURE 6), box 1 10. The jump to box 1 10 marks 

the start of the next or subsequent pair of DR, SEC Functions. The functions of block 
r ^ 15 1 10, path 1 52 and block 150 are identical in operation to those of blocks 90, path 132 

P and block 130 discussed above. The routine then operates through the next traverse 

and again increments the PF index registers. 

The SEC function can operate concurrently with the DR function to process 

the file or set of files most recently generated by the DR function. The SEC function 
20 will for example reprocess the values from the last survey interval files while the DR 

function records values during the current survey interval. 

While the invention has been explained in connection with several 

embodiments, it is intended that the appended claims be used to limit the scope of this 

patent. 



